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Abstract 

The consistent theory of the Heisenberg quantum antiferromagnet in the disordered phase 
with short range antiferromagnetic order was developed on the basis of the path integral for 
the spin coherent states. We have presented the Lagrangian of the theory in a form which is 
explicitly invariant under rotations and have found natural variables in the term of which one 
can construct a natural perturbation theory. The short wave spin fluctuations are similar to 
the spin wave theory ones, and the long wave spin fluctuations are governed by the nonlinear 
sigma model. We have also demonstrated that the short wave spin fluctuations have to be 
considered accurately in the framework of the discrete version in time of the path integral. In 
the framework of our approach we have obtained the response function for the spin fluctuations 
for the whole region of the frequency u and the wave vector k and have calculated the free 
energy of the system. 
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The theory of the two-dimensional Heisenberg antiferromagnet (AF) has attracted great 
interest during the last years in connection with the problem of AF fluctuations in copper 
oxides 0, ||. The approach of these papers was based on the sigma model, which de- 
scribes the long wave fluctuations of the Heisenberg AF in the paramagnetic phase with a short 
range antiferromagnetic order. The sigma model is the continuum model for the unit vector 
n(i, r), n 2 = 1 in the 1 + 2 time and space dimensions |], [/J. As a long wave theory, the sigma 
model can make a lot of physical predictions such as the structure of the long wave fluctuations 
and the magnitude of the correlation length 0, ||, |J . But up to now a consistent theory of the 
spin fluctuations for the quantum AF (QAF) with short range AF order was absent. This is 
just the topic of this paper. 

Our approach to the description of the QAF is based on the functional integral for the 
generalized partition function (GPF) in terms of spin coherent states. We introduce the concept 
of invariant spin coherent states and on this basis we formulate the theory. 

We define the invariant spin coherent states (SCS) with the help of relation: 

|n;m >= exp(— i<pS z ) exp(— i9S y ). exp(— iif)S z )\ss > . (1) 

Here, the state \ss > is the state of spin s with the maximal spin projection s. The unit 
vectors n and m are orthogonal: n 2 = 1, m 2 = 1, n ■ m = 0. 9,(p are the Euler angles of 
the unit vector n = (cosy? sin 9, sin (p sin 9, cos 9). The dependence on the vector m is included 
in the angle ip only, which, in fact determines only the phase factor in the SCS (jl|). We 
can choose the angle ip in some special manner which distiguishes this definition from the 
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standard one ||: tanip = —k z /m z , where the vector k = [n x m]. This choice has a clear 
geometrical interpretation. The transformation (|IJ) rotates the reference coherent state which 
is characterized by the vectors n = (0,0, 1) and m = (1,0,0), into the SCS (H). ^From this 
geometric interpretation it is obvious that upon changing the vector n into the vector n' with 
the help of some rotation a, we have, 

| an; am >= U(a) |n; m > (2) 

without the phase factor which was introduced and discussed by Perelomov ||. It seems that 
the vector m is an artificial one. However, for the problem of the QAF it has some real meaning. 
We consider the spin system which is described by the Heisenberg Hamiltonian: 

H H e l = J - E S, - S,/, § r S, = s(s + l), (3) 

1 i,i'=<i> 

where Si are the spin operators; the index I runs over a two-dimensional square lattice; the 
index I' runs over the nearest neighbors of the site I; J > is the exchange constant which, 
since it is positive, corresponds to the AF spin interaction; and s is the magnitude of spin. The 
most efficient method of dealing with a spin system is based on the representation of the GPF 
Z or the generating functional of the spin Green functions in the form of a functional integral 
over spin coherent states 

Z = Tr [exp (~PH)] , (3 = 1/T, (4) 

/oo roo 
•••/ £>/i(n a ,n b )exp(A(n a ,n 6 )), (5) 
-oo J— oo 

D^(n a ,n b )= J] ^±l5(nl(r,l)-l)dn p (r,l) (6) 

p=a,b;r,l 

where T is the temperature, r is the imaginary time, and A(n) is the action of the system. In 
the continuum approximation, which is valid in the leading order in l/2s the expression of the 
action A(n) is simplified 

A(n a , n 6 ) = - / V C tot (r, l)dr, C tot (r, I) = C kin (r, I) + H{t, I), (7) 
Jo | 

d 

C kin (T,l) = E < n P (T,l);m p (T,l)\-^-\n p (T,l);m p (T,l) >, (8) 

p=a,b ° T 

H(r,l) = Js 2 E n a (r,l)-n b (r,l'). 

i'=<i> 

The idea of the short range AF order was used in Eqs. (|]-|), and we split our square lattice into 
two AF sublattices a and b. For the kinetic part of the action (which is highly nonlinear) 
we use the concept of invariant coherent state parametrized by arbitrary vectors m a) &. 

In our case we can define these vectors in the following manner m a) &: m a = -R a b[n^ — (n^ ■ 
n a )n a ], m b = R a b{n a — ( n a • n b) n b], and R ab = [1 — (n a ■ nf,) 2 ] -1 / 2 , and the invariant coherent 
states have a clear meaning. Substituting these expressions for m ab into Eq. (g) for C kin we 
have also the invariant expression C kin 

C-kin = Z ( n ar/ _ n fcr/) ' [ n ar« x n br«] • (9) 

t — n aT ; ■ Ilhri 
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Now we can introduce new more convenient variables fl(r, I) and M(r, /) which realize the 
stereographic mapping of a sphere: 



±n (1 - M 2 /4) - fn x Ml „ 
n a , b = 1 i + M 2/4 ' "- M = °- ( 10 ) 

In terms of these variables the total Lagrangian Cum — £un + W has the final form 

£ *» = TX^TZ' H = J, 2 E {^-^[(l-M 2 /4)(l-M'V4) (11) 

1 + 1V1 /4 «'=</> 

-M • M'] + n • M' O' • M}(1 + M 2 /4) _1 (l + M /2 /4)~\ 

where O = O r ;, Q' = f2 T //, M = M T ;, M' = M T j/. After this change of variables the measure 
of integration D/j,(n) (|10D becomes 

^(n) = n {2 7J ) \l+mylr 6 (" 2 " 5 • M ) ( 12 ) 

where the product in ([12]) is performed over the AF (doubled) lattice cells. 

The variable f2 is responsible for the AF fluctuations and the variable M for the ferro- 
magnetic ones. The ferromagnetic fluctuations are small according to the parameter l/2s and 
therefore one can expand the Lagrangian Cqm (|i~l"D over M. The vector of the ferromagnetic 
fluctuations M plays the role (up to the factor 2s) of the canonical momentum conjugate to 
the canonical coordinate f2. The term of first order in M coincides (after change of variables) 
with previous results pi [|. 

From Eq. ([[]) one can easily extract the quadratic part of the total lagrangian in the variables 
fi and M, C quad , 

C qU ad = 2is(M ■Cl) + Js 2 J2 - n ■ n ' + m2 + M ■ M '] > ( 13 ) 

i'e<i> 

The Lagrangian C gua d (|13D is very simple but the measure D[i (0) is not simple due to the 
presence of two delta- functions. Therefore we cannot simply perform the Gaussian integration 
over the fields f2 and M. To solve this problem we shall use the method of the Lagrange 
multiplier A together with the saddle point approximation |6|, [7]] to eliminate S(fl 2 — 1). As a 
result, we shall have an additional integration over A with the additional Lagrangian 

C x {r,l) = [iA(r,0 +fi 2 /2J][n 2 (r,l) - 1], (14) 

where fj,Q is the primary mass of the Q field, and J = Jsz. 

To eliminate 5 (CI ■ M) we shall use some kind of Faddev-Popov trick 0. As a result of this 
trick: (1) the factor 5 (ft ■ M) disappears from the measure (p)); (2) M M tr = M-fi(fJ-M) 
in the Lagrangian ([11]); (3) an additional contribution to the action appears, the Lagrangian of 
which Cgaug is conveniently chosen in the form 

£ 9 aug = Js 2 £ [(ft-M) 2 + (ft-M)(ft'-M')], (15) 

z'e</> 
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such choice kills the major dependence on fl in the Lagrangian ([13]) which appears due to 
substitution M =>- M tr ; (4) in the measure of the integration in (0) the additional factor 
(det(B gaug )) 1 ^ 2 , where the operator B gaug is just the operator in the quadratic form in the 
variable (ft ■ M) in ([151) . In this way, the expression ( |I3|) for L qua d is valid in the leading order 
with respect to l/2s. The final expression for the total Lagrangian is L to t = £>q.m + £-gau g + £a 
(0,0,0). 

Now, from the quadratic part of the total Lagrangian Ltot one can find the Green functions 
of the ft and M fields 

L q = u 2 + ^ 2 k , u& = P k Q k = (1 - 7k)^ 2 + (1 + 7k)^/2- 
(Qk, fk) = J(l ± 7k), 7k = (1/2) (cosfea) + cos( V)). 

/FromEq. (ff6|) one can calculate the parameter of the spin wave nonlinearity of the theory: 
< M 2 tr >= (l/2s)C M (T), where C M (T) = 0.65075 for T < J, and C M (T) = 1.48491T/J for 
T > J. 

We also have the saddle point condition for the A field < ft 2 >= 1 which is the most 
important constraint of the theory which determines its phase state: 

i=<«*>=iv£G? = ^ Y. E^ct- (17) 

q ZS u,=2nnT k U + ^Ok 

l=<tt 2 >=f £^(l + 2n 0k ), (18) 
2s k 2u; k 

where n k = (exp(co»ok/T) — l)^ 1 is the Plank function. The right hand side of Eq. (0) 
contains two terms. The first term Qk/2o;ok is responsible for the quantum fluctuations of the 
fl fields. The second term Qk^ok/^ok is responsible for the classical thermal fluctuations of the 
fl fields. The role of these two terms is quite different. The quantum fluctuations are small 
according to the parameter of perturbation theory l/2s and, for the basic approximation, they 
can be neglected. The thermal fluctuations can be considered in the continuum approximation 
which leads to the well known Jl|, |2|, H zero order expression for fiQ 

Ho = Texp ^~ 2 y^ 1 , € = nc a /iJi, (19) 

where £ is the correlation length. From Eq. ([19]) the important conclusion is followed: in the 
regime of the weak coupling the correlation length £ is much larger than the lattice constant a. 
To close the theory it is helpful to define the polarization operator 11(g) of the fl field 

A A?uad = -~£A*(g)II(g)A(g), (20) 

and the Green function of the A field is 11(g) _1 . In the lowest approximation 11(g) is simply a 
loop from two Green function G 

U ( q ) = 2NTj2G n (q')G n (q-q'). (21) 
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The main contribution in l/2s in (|2T| ) is from the thermal fluctuations even at low temperatures 
T, because the integral strength of such fluctuations is fixed by the saddle point condition (|17|) 
and does not depend on the temperature. The explicit form for Tlo(q) one can get in two 
limiting cases hq ^> T and hq -C T, where q 2 = uj 2 + c 2 k 2 . In the first case the momentum 
q' ~ T/c s <C q, and we can separate summation and integration over q' and put q' = in 
G (q — q') in (0). The result is extremely simple 

n (g)=4^(g)= 2 f ( 2 1+7 2 k) g»fc T , k T = T/c s . (22) 

Notice, that it exceeds the quantum contribution in ([24]) Ilo(g) = iV/4g by the large parameter 
WsJ'/Nq. For small g <C c s /a and q <^kr our results coincide with 0. 

The dynamical spin susceptibility ^ (cu, k) for all values of uj and k can be calculated. In 
the lowest order over l/2s we can use the lowest order relation 

n(n(r, 0, M(r, /), r, Z) ~ e ial ^n(r, - [n(r, /) X M(r, /)], (23) 

where q^F — (7r/o, 7r/a) is the AF vector (|T0|). Calculating the average of two vectors n from 
(HI) we get the dynamical spin susceptibility as a sum of two terms Xij{ u i k) = ^ijlXA^, k) + 
XF(^,k)]. The spin susceptibility XA(u,k) is responsible for the AF fluctuations. It is propor- 
tional to the Green function analytically continued and shifted by the AF vector qaf 

i \\ Js 2 z(l + ^) 

Xa(w, k = — 777 2 2 T • K > 24 

2(u;^ — cj>Q k » + «u;o) 

where k* = k — q^i?. For the ferromagnetic spin susceptibility Xf(^ ; k) we have a loop expres- 
sion which can be calculated on the basis of the thermal fluctuation domination, as a result we 
have for q>kx 

X> w, k ^ — irrG (q) = - 1 2 ; kJ 25 

The theory of the spin fluctuations in the disordered QAF at sufficiently low temperature 
T J allows for the scale separation. In this case k^ n/a and the thermal fluctuations 
can be considered by the "renormalized classical" manner @. The magnitude of the quantum 
fluctuations at q < kx is small in comparison with the classical fluctuations. In this situation 
the parameters of the effective long wave, low frequency sigma model are renormalized by 
the quantum fluctuations . This renormalization is performed with respect to the parameter 
l/2s, but the interaction of the thermal fluctuations with the scales |k| < kx and uj < T is 
over parameter 1/N, where iV is the number of components of the n field of the long wave, 
low frequency nonlinear sigma model. This picture follows directly from the approach of this 
paper. 

Unfortunately, the continuum approximation in time is not working when we are calculating 
corrections to the basic approximation. The reason for last observation lies in the canonical 
structure of the Lagrangian ( |TT| ) and the Green function flOp: the sums over uj including this 
Green function are ambiguous and must be defined at the final time step A. We only describe 
the basic approximation and results shortly. They can be obtained on the basis of the SCS (|T|). 
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Instead of the expression (0) for the action A(n) we shall use more accurate expression 

A (n) = - E E A [A;n(j, + H(j, I)] , (26) 

j=0 I 

where r = jA, and AN T = (3. Here, £kin(j, I) consists of two parts Ckin = £mod + £ P ha- The 
first term is pure real the second term is pure imaginary. 



AC mod = -sin [(1 + na- n n )(l + n 6 • n 6 )/4] 



(27) 



where n p = n p (j, l),n p = n p (j + 1, /), for p = a,b. 
The Lagrangian C p h a is not so simple 



ACpha 



E J" 



p=a,( 



(28) 



where the quantity R p = n p -m p +in p -k p , where vectors n, m, k were defined at the introduction 
of the SCS. 

Expanding £ p h a in the vector M has a rather complicated form but one can prove that it 
is regular and contains only odd powers of M. 

The Hamiltonian 7~C(n) can be obtained on the basis of the following relation 



<S(n, n) 



<n|S|n> n + n — z[n x n] 



< n n > 



n ■ n 



(29) 



for the matrix element of the spin operator S. If we substitute them into the matrix element 
of the Heisenberg Hamiltonian we obtain 



W(n) = Js 2 E S(n,n)-S(n',n'), 

i'e<i> 



(30) 



It is assumed that all vectors n p , m a , k p for p = a,b entering in Eqs. ( p7| - |30|) are functions 
of the dynamical variables O and M according to Eq. fllT^) . 

By expanding the Lagrangians £ m od (P?D , C p h a (^8|), and the Hamiltonian (130|) in the vector 
M up to second order we get 

AC kin = s[l - Q, ■ n + M 2 - M ■ M + i(Q • M - ft ■ M)}, (31) 

n = Js 2 E [Q • n - n • n' +m • m + m • m' - i(o • m - n • m)]. 

i'e<i> 

Now we present the result for the free energy Fqaf of QAF which was obtained in the lowest 
order in l/2s on the basis the Lagrangian (|3~T|). Fqaf = FnM + Fx, and 

(32) 



Fqm = ~N S J + 2N S E {^ok/2 + T ln[l - exp(-aWT)]} 

k 



A 



E 1 * 



o.'k 



g(o; 2 + ^ 2 k )n (g) 
2J(l + 7 k) 
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Here Fqm is the free energy of the f2 and M fields, F\ is the free energy of the A field, 



2N S is the number of the lattice sites, and the polarization operator n (q r ) is defined in (|2T|). 
The temperature dependent part of the free energy (j32|) at small temperatures T <C J is 
proportional to Faf ~ NgT 3 /^. Such contribution has two origins: one from Fqm and another 
one from F\. 

In the framework of the discretized time approach, the first order corrections to the Green 
functions of the f2 and M fields were calculated and the expression for the correlation length 
in the next order in l/2s was obtained which lies completely in the framework of the concept 
of scales separation. This results will be discussed in detail in a more complete publication. 
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